Optimal spin squeezing inequalities detect bound entanglement in spin models 
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We determine the complete set of generalized spin squeezing inequalities. These are entanglement 
criteria that can be used for the experimental detection of entanglement in a system of spin-i 
particles in which the spins cannot be individually addressed. They can also be used to show the 
presence of bound entanglement in the thermal states of several spin models. 
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Entanglement lies at the heart of many problems in 
quantum mechanics and has attracted increasing atten- 
tion in recent years. However, in spite of intensive re- 
search, many of its intriguing properties are not fully un- 
derstood. For example, it has been shown that there are 
entangled states, from which the entanglement cannot 
be distilled again into the pure state form, even if many 
copies of the state are available jlj. The existence of 
these so-called bound entangled states has wide-ranging 
consequences for quantum cryptography Q and classi- 
cal information theory [3j. Since entangled states that 
are not recognized by the separability criterion of the 
positivity of the partial transpose (PPT) [3| are bound 
entangled, such states also serve as a test bed for new 
separability criteria [5|, [6|, |7fl. However, bound entangled 
states are often considered to be rare, in the sense that 
they do not occur under natural conditions. 

In physical systems such as ensembles of cold atoms [8[ 
or trapped ions Q, spin squeezing 1^, U| is one of the 
most successful approaches for creating large scale quan- 
tum entanglement. Since the variance of a spin compo- 
nent is small, spin squeezed states can be used for re- 
ducing spectroscopic noise or to improve the accuracy 
of atomic clocks [l^, Hl||. Moreover, if an 7V-qubit state 
violates the inequality [121 ] 




FIG. 1: The polytope of separable states corresponding to 
Eqs. ([2} for N = 6 and J = 0. S corresponds to a many body 
singlet state. 



individually addressed, and only collective quantities can 
be measured. In Ref. [l3| a generalized spin squeezing 
criterion was presented detect ing the presence of two- 
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where J; := | Efc=i a i f° r ^ = x >Vt z are the collective 
angular momentum components and er, are Pauli ma- 
trices, then the state is entangled (i.e., not separable), 
which is necessary for using it in quantum information 
processing applications 

Recently, several generalized spin squeezing criteria for 
the detection of entanglement appeared in the litera- 
ture [HI, 0, EH and have been used experimentally . 
These criteria have a large practical importance since in 
many quantum control experiments the spins cannot be 



qubit entanglement. In Refs. [14 . |15| other criteria can 
be found that detect entanglement close to spin singlets 
and symmetric Dicke states, respectively. These entan- 
glement conditions were obtained using very different ap- 
proaches. At this point two main questions arise: (i) Is 
there a systematic way of finding all such inequalities? 
Clearly, finding such optimal entanglement conditions is 
a hard task since one can expect that they contain com- 
plicated nonlinearities. (ii) How strong are spin squeezing 
criteria? Can they detect entangled states that are not 
detected by the PPT criterion or other known entangle- 
ment criteria? 

The goal of this Letter is twofold. First, we give a 
complete set of spin squeezing inequalities based on the 
first and second moments of collective observables. Sec- 
ond, we use them to show the presence of multipartite 
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bound entanglement in several spin models in thermal 
equilibrium. In particular, we consider bound entangle- 
ment that has a positive partial transpose with respect 
to all bipartitions. 

We can directly formulate our first main result: 
Observation 1. Let us assume that for a physical 
system the values of J := (( J x ), ( J y ), ( J z )) and K :— 
((J x ), (Jy), (Jz)) are known. Violation of any of the fol- 
lowing inequalities implies entanglement: 

(Jx) + (Jy) + (Jz) < N(N + 2)/±, (2a) 

(AJ X ) 2 + (AJy) 2 + (AJ Z ) 2 > N/2, (2b) 

(J?) + {J]) - N/2 < (N-l)(AJ k ) 2 , (2c) 

(N-VliAjtf + iAJj) 2 ] > (J 2 )+N(N-2)/4,(2d) 

where i,j,k take all the possible permutations ofx,y,z. 
The proof is given in the Appendix. 

For any value of J these eight inequalities define a poly- 
tope in the three-dimensional (( J 2 ), ( J 2 ), ( J 2 })-space. 
Observation 1 states that separable states lie inside this 
polytope. For the case J = and N — 6 the polytope is 
depicted in Fig. [TJ Such a polytope is completely char- 
acterized by its extreme points. Direct calculation shows 
that they are 
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where k := (JV - 1)/N. The points A y/ 
obtained in an analogous way. 

One might ask whether all points inside the polytope 
correspond to separable states. This would imply that 
the criteria of Observation 1 are complete, that is, if the 
inequalities are satisfied, then the first and second mo- 
ments of Jfc do not suffice to prove entanglement. In 
other words, it is not possible to find criteria detecting 
more entangled states based on these moments. Due to 
the convexity of the set of separable states, it is enough 
to investigate the extreme points: 

Observation 2. For any value of J there are separable 
states corresponding to Ak ■ For certain values of J and 
N there are separable states corresponding to points Bk- 
However, there are always separable states corresponding 
to points B' k such that their distance from Bk is smaller 
than 1/4. In the limit N — > oo for a fixed normalized 
angular momentum j := J /(N/2) the difference between 
the volume of polytope of Eqs. 0) and the volume of set 
of points corresponding to separable states decreases with 
N at least as AV/V cx N~ 2 , hence in the macroscopic 
limit the characterization is complete. 

Proof. A separable state corresponding to A x is 



Here |V>+/-) are the sm gl e qubit states with Bloch vec- 
tor coordinates ((a x ), (a y ), (a z )) = (±c x , (J y )/J, (J*)/ J) 
where J := N/2 and c x := yjl - ((J y ) 2 + (J z ) 2 )/J 2 . 
The mixing ratio is defined as p := [1 + (J x ) / (Jc x )]/2. 
If M := Np is an integer, we can also define the state 
corresponding to the point B x as 
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If M is not an integer, we can approximate B x by tak- 
ing m := M — e as the largest integer smaller than M, 
defining pf := (l-e)(|^ + )(i/' + |)® m ®(|V'-}(V'-l) 8(Ar_m) + 
e(|i/' + }(V'+l) 8(m+1) ®(IV'-}(V'-l)® (Ar_m_1) -This state has 
the same coordinates as B x , except for the value of (J 2 ), 
where the difference is c 2 (e — e 2 ) < 1/4. The depen- 
dence of AV/V on N can be studied by considering the 
polytopes in the ((J x ), (Jy), ( J z })-space corresponding to 
(Jk) = jk x N/2, where jk are the normalized angular 
momentum cordinates. As N increases, the distance of 
the points Ak to Bk scales as N 2 , hence the volume of 
the polytope increases as ./V 6 . The difference between the 
polytope and the points corresponding to separable states 
scales like the surface of the polytope, hence as iV 4 . □ 

Now we consider already known entanglement criteria 
and show how they can be derived from our theory. This 
can be done by showing that for any J the points Ak and 
Bk satisfy them. 

Case 1. The standard spin-squeezing inequality is 
Eq. |T]) from Ref. 12]. This inequality is valid for all 
Ak and Bk, for B x even equality holds. 

Case 2. For separable states (J 2 ) + (J 2 ) < (N 2 + 
N)/A holds [IB}, as can be proved in the same way. This 
can be used to detect entanglement close to the iV-qubit 
symmetric Dicke states with N/2 excitations. 

Case 3. Separable states fulfill Eq. (|2b[) which has 
already been shown in Ref. [3] • It is maximally violated 
by a many-body singlet, e.g., the ground state of an anti- 
ferromagnetic Heisenberg chain. 

Case 4- For symmetric states it is known that (J 2 ) + 
(Jy) + (Jz) = N ( N + 2)/4 [13]. From this and Eq. jScJ 
one can directly derive 1 - 4(Ji) 2 /N 2 < 4(A Ji) 2 /N from 
Ref. [H . 

Next, it is interesting to ask what kind of entanglement 
is detected by our criteria knowing that they contain 
only two-body correlation terms of the from (a[ <Jju ) 
and do not depend on higher order correlations. In 
fact, all quantities in our inequalities can be evaluated 
based on knowing the average two-qubit density matrix 
Pav2 := jv(jv-i) ^fijtj Pir Do our cr i teria simply detect 
entanglement of the two-qubit reduced state of the den- 
sity matrix? We will now show that the criteria Eqs. ([2]) 
can detect entangled states that have a separable two- 
qubit density matrix. Even more surprisingly, they can 
detect bound entanglement in spin systems. While in the 
following we will use Eqs. Q for the theoretical analy- 
sis of spin models, we stress that Eqs. ([2) can also be 
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PPT[1|234] -— 
PPT[12|34] — 
PPT[13 24] — 
Permutations — 
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FIG. 2: Entanglement properties of the spin model with the 
Hamiltonian Eq. ([5]). The critical temperatures for several 
entanglement conditions are shown as a function of the next- 
to- nearest neighbor coupling J2. For details see text. 



used for the experimental detection of entanglement in a 
realization of these models in physical systems in which 
the collective angular momentum can be measured (e.g., 
Ref. 0). 

Let us first consider four spin-1/2 particles, interacting 
via the Heisenberg-type Hamiltonian [l7T | 



h = &k&k+i + ^2(^1 <? 3 + CT2CT4) 



(•5) 



k=i 



where a = (a x , cr yi a z ). For the above Hamiltonian, we 
compute the thermal state g(T, J 2 ) oc exp(— H / kT) and 
investigate its separability properties. Hamiltonians of 
the type Eq. are by no means artificial: They are 
used to describe cuprate and polyoxovanadate clusters 



17, 18]. For several separability criteria we calculate the 
maximal temperature, below which the criteria find the 
states entangled. The results are summarized in Fig. [2] 
For J2 ^ —0.5, the spin squeezing inequality Eq. (|2b]) 
is the strongest criterion for separability. It allows to 
prove the presence of entanglement even if the state is 
PPT with respect to all bipartitions This implies 
that the state is multipartite bound entangled: No pure 
entangled state can be distilled from it [l9(. Note that 
introducing the next-to-nearest neighbor coupling made 
the PPT entangled temperature range larger. 

For comparison, we investigated the computable cross 
norm or realignment criteria (CCNR, ||) corresponding 
to all bipartitions, all the other criteria based on permu- 
tations [201 ] , and the criterion based on covariance matri- 
ces ■ None of them is able to find PPT entanglement in 
our spin system. Finally, we studied for each bipartition 
the separability test of symmetric extensions [2 11 ] that 
is strictly stronger than the PPT criterion. The critical 
temperatures, however, coincide within numerical accu- 
racy with the ones from the PPT criterion, giving strong 
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2.56 


3.46 
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3.34 


3.19 


3.32 


3.24 



TABLE I: Critical temperatures for the PPT criterion and 
Eqs. ([2]) for Heisenberg and XY spin chains of various size. 



evidence that g is indeed separable for the bipartitions. 
Indeed, we will see later that in some spin models, the 
spin squeezing inequalities signal the presence of entan- 
glement even for states that are separable with respect 
to all bipartitions. 

After small spin clusters, we consider larger spin sys- 
tems. Using Eqs. ([2]), we find bound entanglement that 
is PPT with respect to all bipartitions in Heisenberg and 
XY chains with a periodic boundary condition with up to 
9 qubits. The critical temperatures are shown in Table 
I. Eqs. ([2]) also detect bound entanglement in Heisen- 
berg and XY models with a complete graph topology 
[22| . Latter is a special case of the Lipkin-Meshkov-Glick 
model [233. In all these cases there is a considerable 
temperature range for which the thermal state is PPT 
with respect to all partitions but not yet separable [24j ]. 
Interestingly, since in the three-qubit Heisenberg model 
the thermal state is invariant under multilateral unitary 
transformations of the type U ® U ® U, for such states 
the PPT condition implies biseparability 25| . Thus, the 
spin-squeezing inequalities can detect bound entangle- 
ment for which all bipartitions are separable. 

Note that the bound entanglement that is PPT with 
respect to all bipartitions is perhaps the most intrigu- 
ing type and the most challenging to detect. No pure 
state entanglement can be distilled from it with local op- 
erations and classical communication, even if arbitrary 
number of parties join. However, an entangled state that 
is PPT with respect to only a single partition is already 
bound entangled since no GHZ state can be distilled from 
it [19( . Such entanglement can be found by the PPT cri- 
terion with respect to a different partition. It is expected 
to appear in many systems since as the temperature in- 
creases, not all the bipartitions become PPT at the same 
temperature [2f|. 

Our study of the spin models has two general conse- 
quences. First, we realize that examination of spin mod- 
els via the partial transposition or the investigation of 
bipartitions does not lead to a full understanding of the 
entanglement properties of condensed matter systems. 
Second, we note that the spin clusters and spin chains 
we studied are models of existing physical systems. Thus 
multipartite bound entanglement that is PPT with re- 
spect to all partitions is not a rare phenomenon in na- 
ture. 



Moreover, based on Ref. [27[, it is possible to connect 
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the variances of collective angular momenta to important 
thermodynamical quantities giving our inequalities a new 
physical interpretation. Let us consider a system with a 
Hamiltonian H and an additional magnetic interaction 



B k J k , where B is the magnetic field. 



Hi ■= E fe = . 
Moreover, assume that H commutes with J x / y / z - Then 
the magnetic susceptibilities are \i '■— (9(Ji)/dBi)\g =Q 
for I — x,y,z and the variances can be written as 
(AJ/) 2 = kTxi- Thus our inequalities can be expressed 
with susceptibilities [28j |. 

Finally, we discuss some further features of our spin 
squeezing inequalities. One can ask what happens, if not 
only (J k ) and ( Jj 2 ) for k = x,y, z are known, but ( Jj) and 
(J 2 ) in arbitrary directions i. We will now show how to 
find the optimal directions x' , y' , z' to evaluate Observa- 
tion 1. Knowledge of (Jj) and (J 2 ) in arbitrary directions 
is equivalent to the knowledge of the vector J, the corre- 
lation matrix C and the covariance matrix 7, defined as 
0,111 Cm := (J k J t + JJ fe )/2 and lkl := C kl - (J k )(Ji) 
for k,l — x, y, z. When changing the coordinate system 
to x',y',z', vector J and the matrices C and 7 trans- 
form as J n OJ, C 1— » OCO T and 7 1— » OjO T where 

is an orthogonal 3 x 3-matrix. Looking at the in- 
equalities of Observation 1 one finds that the first two 
inequalities are invariant under a change of the coordi- 
nate system. Concerning Eq. (j2c|) . we can reformulate it 
as (J 2 ) + (J?) + (4) - N/2 < (JV - l)(AJ fe ) 2 + (J 2 ). 
Then, the left hand side is again invariant under rota- 
tions, and we find a violation of Eq. ()2c|) in some di- 
rection if the minimal eigenvalue of X := (JV — 1)7 + C 
is smaller than Tr(C) — N/2. Similarly, we find a vio- 
lation of Eq. (|2d|) if the largest eigenvalue of X exceeds 
(JV - l)Tr(7) - N(N - 2)/4. Thus the orthogonal trans- 
formation that diagonalizes X delivers the optimal mea- 
surement directions x',y',z' .'SO . 

In summary, we presented a family of entanglement 
criteria that detect any entangled state that can be de- 
tected based on the first and second moments of collective 
angular momenta. We applied our findings to examples 
of spin models, showing the presence of bound entangle- 
ment in these models. 
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the FWF, and the Spanish MEC (Ramon y Cajal Pro- 
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Appendix — Proof of Observation 1. Fully separable 
states are of the form p = YliPlPi ® Pi ® •■■ ® /°j > 
where Ei-P' = 1 an< 4 Pi > 0. The variance, defined 
as (A A) 2 := (A 2 ) — (A) 2 , is concave in the state thus 
it suffices to prove that the inequalities of Observation 

1 are satisfied by pure product states. Based on the 
theory of angular momentum, Eq. (|2ap is valid for all 



quantum states. For Eq. (|2b[) one first needs that for 

product states (AJ fc ) 2 = JV/4 - (1/4) EjH^) 2 holds > 
then the statement follows form the normalization of the 
Bloch vector. Concerning Eq. l[2c|) . we have to show 
that 2) := (JV — 1)(AJ X ) 2 + N/2 - (J 2 ) - (J 2 ) > 0. 

Using the abbreviation x, — {o~x^),yi = (o~y ), etc. this 
can be written as 2) = (JV - l)[JV/4 - (1/4) E> ^f] ~ 
(1/4) (WW + = (JV - l)[JV/4 - (1/4) Y,i - 
(l/4)[(Ei y,) 2 + (Ei ^) 2 ] + (1/4) E,(2/ 2 + zl)- U^ng the 
fact that (Ei s i) 2 < ^Ei s f' an< 4 l ne normalization of 
the Bloch vector, it follows that 2} > 0. Eq. (|2d[) can be 
proved in the same way. □ 
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